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I. INTRODUCTION 

The calculation of weak matrix elements, such as the fs decay constant, is one of the most important applications 
of lattice QCD, because it is this nonperturbative physics that obscures possible physics beyond the standard model. 
Although lattice QCD offers a nonperturbative method of calculating weak matrix elements from first principles, in 
practice a perturbative calculation is also required to extract the continuum number. Lattice perturbation theory 
is much harder than continuum perturbation theory, because there is less symmetry on the lattice that can be used 
to simplify expressions. This means that lattice perturbative quantities are typically only known to one loop in 
perturbation theory. It seems important to extract as much information from the one-loop calculation as possible. 

In the past, one-loop lattice perturbation theory was considered unreliable, because for many quantities the one- 
loop corrections are large if one uses the bare lattice coupling constant, and the agreement with results that could be 
calculated numerically was the very poor. Lepage and Mackenzie [0 then showed that the use of a more “physical” 
coupling, calculated at a judiciously chosen scale (generically referred to as q*) and the resummation of gluon-tadpole 
contributions to the perturbative expansion improve the situation for many quantities. 

In this paper we calculate the scale q * using the Lepage and Mackenzie prescription for the local axial and vector 
currents (as well as the mass times the pseudoscalar density), with and without tadpole improvement. Some of our 
results were also obtained in Ref. Q; however that paper did not study tadpole-improved perturbation theory. In 
addition, the current work shows how to deal with certain ambiguous (though numerically small) contributions to the 
calculation in Ref. (Q. 

There are a number of nonperturbative renormalization techniques that have been used in numerical computa¬ 
tions [|] |5). These techniques are very interesting, but they have not as yet been widely used in large scale simulations. 
The calculation of the nonperturbative renormalization constants required in the calculation of fs in the simulations 
by the MILC collaboration |j|, for example, would require, at the least, a substantial investment in computer time. 
There are also conceptual problems in generalizing existing nonperturbative methods to the case of large quark mass. 
Furthermore, it is an important consistency check on the numerical calculation of renormalization factors, that they 
should agree with the weak-coupling results at small enough coupling, because there are potential systematic errors, 
such as lattice-spacing artifacts, in the numerical results. This motivates trying to extract as much information as 
possible from lattice perturbation theory calculations. 

The outline of this paper is as follows: in section 2 we calculate the renormalization constants from Ward identities, 
while in section 3, they are calculated by matching on-shell matrix elements to the continuum. In section 4, these 
renormalization constants are tadpole improved, and in section 5, we calculate the scales q * for each of them with and 
without tadpole improvement. In section 6, we discuss the stability of our results under changes in the scale-setting 
prescription. We end by providing some insight into the fact why tadpole improvement does not lower the scale q * 
for these quantities as much as one might naively expect. 
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II. Z A , Z mP AND FROM WARD-IDENTITIES 


First, we will derive the renormalization factors from the Ward identities |t],|]]. We consider the following Ward 
identity (WI) for the nonsinglet flavor axial current A M = A“T a ( T a are the hermitian, traceless flavor generators; we 
work in euclidean space): 

(52 d pK( x )^ z )q{y)) = 2 m{P a (x)q(z)q{y)) (2.1) 

A 4 

-{q(z)q(y))T a ^ 5 S(x -y)- T a j 5 (q(z)q{y))S{x - z) , 


with 


A pi x ) = q(xh/i.l5T a q(x) , 

P(x) = q(x)j 5 T a q(x) . 


( 2 . 2 ) 

(2.3) 


This WI holds for any regulator which respects the (nonsinglet) axial flavor symmetry, such as naive dimensional 
regularization (NDR), in which y 5 anticommutes with 7 ^ for all y = 1,... ,d. Since is the only dimension-three 
axial vector operator that can be built out of quark and gluon fields, it does not mix with other operators. Taking 
m = 0 in eq. & we see that does not require renormalization (the RHS of eq. is finite after renormalizing 
the quark fields q(z) and q(y)). Thus, if Za is defined to be the renormalization constant for which preserves the 
WI (2T), we have Za = 1 in any scheme, such as NDR, which preserves the axial flavor symmetry. The operators 
mP a also do not mix with any other operators, and therefore Z m p (defined to be the renormalization constant of 
the operator mP a ) may be determined by considering the case m ^ 0 after Za has been defined. In NDR, we have 
ZmP — 1 * 

Wilson fermions break the axial flavor symmetry, and no conserved axial currents exist, so that the WI introduced 
above does not hold for any axial current on the lattice. Here, w e are interested in the local («.e., single-site) 
axial current and pseudo-scalar density as defined in eqs. (|2.2 .2.3), and we would like to find the multiplicative 
renormalization constants Za and Z m p such that eq. (|2.1) holds in the scaling region: 


z a( 52 d ^ A ^ x ) ( l(z)q(y)} = 2Z mP (M - M c )(P a (x)q(z)q(y)) 

~(q(z)q(y)}T a ~/ 5 6(x - y) - T a ^ b {q{z)q(y))8(x - z) 


(2.4) 


Here M = m + 4r (r is the Wilson parameter and m is the bare mass) is the coefficient of the single-site term in the 
Wilson-Dirac lagrangian, and M c is its critical value. On the lattice, we may write 


<£ A ^A;(x)q(z)q(y)) = 2 (M - M c )(P a (x)q(z)q(y)) 


(2.5) 


-(q{z)q(y))T a 'y 5 5(x -y)~ T a j 5 (q(z)q(y))8{x - z) 

+(E a (x)q(z)q(y )) , 

where is the backward difference operator A fJ/ f(x) = f(x) — f(x — y), and E a is the “evanescent” operator 

E a {x ) = Y, A„H“(z) - 2(M - M c )P a (x) (2.6) 


in which 


+q{x)”/ 5 T a {D + M)q(x ) + q(x)(D + M)^T a q{x) , 


Dq{x) =~Y [('V - r )Up{x)q{x + y) - (7m + r ) u t( x - y)q{x - y)] 


(2.7) 


(D + M is the Wilson Dirac operator). 

Going to momentum space by multiplying eq. (2.5) by Ylxyz e^pi—ikx — ipz + ip'y), and amputating the external 
quark lines, we find, at one-loop order in g 1 


(E a (k)q(p)q(p f )) r p - l0 ° P = C F 9 2 S(p + k - p')T a I(p,p') 


( 2 . 8 ) 


2 







with Cf = — 2 ^ 1 the quadratic Casimir for SU(N), and 

I(P,P') = x (2-9) 

fl •’Z 

(r~S(p + £) 7 , 75 S(p' + £)T~ e^+r'+V- 
+ T~S(p + £) lvlh S(p' + £)T+ 

+ r +S(p + ih^ 5 S(p' + £)T~ 

+ T+S(p + ^) 7,75 S(p' + £)T+ e-^+v'+V^ 

+2(M - M c ) J2 [ D{£) (r pS(p + £h 5 S(p' + £)T~ e^'+^ 

fj, Je \ 

+ T~S(p + £)j 5 Sip 1 + £)T+ e^ p - p '^ 


+ T+S(p + £) lb S(p' + £)T~ e-^-P'^ 

+ T+S(p + £)^S(p' + £)T+ e - 

-S (p, m )75 - 7 5 S(p 7 , to) , 

£(p,H = E /^(0 e i(2 ^ } " - r+s(p + £)r+ e-^ 2 ^)" 

-r-s(p + £)r+-r+5( P + £)r- 


( 2 - 10 ) 


2 smp M 7 M - - cosp M 


Here E (p,m) is the one-loop fermion self-energy (defined to be equal to the sum of diagrams, with no overall mi- 

r " i -i 

nus sign), and S(p) = i 7 M sinp^ + m + r (1 — cosp^) is the tree-level fermion propagator. The gluon 
propagator in the Feynman gauge is 6 l _ lu D(p) with D~ 1 (p) = 4 V t sin 2 (p M /2). The matrices Tj are defined by 


't = 2 ^ ±r ) ■ 


( 2 . 11 ) 


Furthermore, we abbreviated J £ = lie ]<* ,d‘ i £. 

We now ob serv e that 7(0,0) = — 275 E( 0 , 0 ) provides the counterterm to remove the 1/a divergence in the contact 
terms in eq. ( |2.5| ) (after amputation of external fermion lines). A straightforward calculation then shows that the 
remainder, I(p,p') — 1(0,0), is finite, and we can take the continuum limit. We expand to linear order in p, p' and 
to = M — M c + 0(g 2 ): 

Cf (I(p,p ) - 1(0,0)) = c A i($' -j >)75 - 2 c mP TO75 , (2-12) 

from which we then obtain the renormalization constants 

Za = 1 g ca j Z mP = 1 g c m p , (2.13) 
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with, for N = 3 and r = 1 


ca = 0.133373(2) , c mP = 0.081419(2) 


(2.14) 


At this point, we pause to comment on this derivation. Generally, in calculating the continuum limit of l attice 
integrals such as I(p,p'), special care is required near the origin in momentum space. Near £ = 0 in eq. ( |2.9|) , one 
may approximate the integrand by its covariant form. In this case, with the routing of the external momenta through 
the fermion line, it turns out that the covariant form of the integrand vanishes identically (as explained below), and 
therefore eq. (2.12) is obtained by a straightforward expansion in p , p' and to. This would not be true had we chosen 
the routing of the external momenta through the gluon line. This will turn out to be relevant for the definition of q *, 
and will be discussed further below. 

A completely analogous analysis may be performed in order to find the renormalization constant Zy of the local 
vector current 


V“(x) = q{xA,jT a q{x) . 

For Wilson fermions, this current is not conserved, and we find 


Z v = l- g 2 c v , Cy = 0.174083(3) . 


Equations ([2.14 2.1C) are in agreement with previous computations 


(2.15) 


(2.16) 


III. Z A , Z mP AND Zv BY ON-SHELL MATCHING TO THE CONTINUUM 


A common approach in the literature is to compute the renormalization constants by matching an on-shell lattice 
perturbative computation to the continuum, where the renormalization constants are known. We present this approach 
here because the value of q* depends on the specific form of the integrands. It is therefore important to check that a 
separate derivation of the integrals will produce a consistent result for q *. 

For definiteness, we focus first on the nonsinglet axial current, A The derivations for the vector (V^) and 
pseudoscalar ( mP ) cases are similar, but there are a few important differences, which we note explicitly. Because A^ 
has no anomalous dimension or mixings, we may write its one-loop matrix element between on-shell quark states i 
and / in regularization scheme S as 

(/ I (A;) s | i) = 7m7 5 T“(1 + g 2 C s A {m, A)) , (3.1) 

where to is the quark mass, A is a gluon mass that has been inserted to regularize the infrared divergences, and spinors 
on the external lines are implicit. 

We then define a on-shell renormalized axial current (A“) sheU 


(A;y™ l = z%(m,\)(A;r 

Z A (m, A) = 1 - g 2 C|(m, A) + . 


(3.2) 

(3.3) 


By definition, the on-shell quark matrix element of (A “) she11 is the same in any scheme and equal to 7 ^ 75 T°, as long as 
to and A are the same as in eq. (3.2). Because there is no operator mixing, this then implies the scheme independence 
of all renormalized Green’s functions with a (A “) she11 insertion. However, we may also define a renormalized current 
(A“) WI using the Ward Identity, eq. (|2.4[): 


(A“) WI = Zi ( A;) s , (3.4) 

where the superscript S specifies the scheme to which eq. ( |2.4| ) is applied. Z A is UV finite; it is free of infrared 
divergences since it is defined off-shell. By dimensional analysis, Z A is therefore also independent of the quark mass 
m. Since renormalized Green’s functions with a (A“) WI insertion are also scheme independent, we have, for any two 
schemes S and S', 


yS yS' 

Z s A (vn,\) Z^ (m, A) 


(3.5) 
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We now specialize to the case S = lattice and S' = NDR. As long as the regulator does not violate the symmetry, 
the Lie algebra of conserved charges fixes the on-shell matrix elements of the corresponding currents, guaranteeing 
that a conserved nonabelian flavor current is not renormalized. Therefore 


ZT R (m = 0, A) = 1 


(3.6) 


since the nonsinglet axial current is conserved when m = 0. Similarly, 

l (to, A) = 1 


7NDR7 


(3.7) 


for all to, A. A simple computation shows, however, that ■Z’( ( pr (to, A ) ^ 1 for general to, A, (including the limit A —> 0 


to fixed). Note that the presence of the gluon mass in eqs. (|3.6| , p.7[) does not violate the flavor symmetries, although 
it would violate gauged symmetries. 

It is useful to see explicitly how eqs. (3.7) and ([b(]) arise at one-loop. At this order, the computations are identical 
to those in QED, and amount the statement that Z\ = Z% 121. If one routes the external momentum p thr oug h the 
fermion line, then the cancellation of the wave-function and vertex renormalizations which results in eq. (|3.7| ) can 
easily be shown from the (Euclidean) identity 


d 


dpv (i i> + to) (i + to) 


(-*7^)77 


(i + to) 


Equation (3.8), coupled with the on-shell definition of the wave-function renormalization, 


Zo = l- 


p 11 <9E 
to dp v 


(3.8) 


(3.9) 




immediately gives Z\ = Z 2 . With the chosen momentum routing, the equality holds at the level of the Feynman 
integrands. On the other hand, if the external momentum is routed through the gluon line, there is no obvious relation 
between the integrands for Z\ and Zi- The equality then appears only after the momentum and Feynman parameter 
integrals are performed and requires a regulator which permits shifts of the loop momentum without the generation of 
boundary terms. This is the reason that the “covariant part” of the lattice integral in Section [0] vanishes only when 
the external momentum is routed through the fermion line. Although complete lattice integrals are invariant under 
shifts, the “covariant part” is not because it is the integral of a continuunr-like integrand over a sphere of radius S 
(the “inner region” of Ref. [[[ 3 ) ) or of radius n (the integration region chosen for the continuum-like integrals added 
and subtracted in Ref. [lj||). The boundary term generated by a shift in the loop momentum thus does not vanish in 
the continuum limit for a linearly divergent integral such as the self-energy. 

For mP , one can relate the mass derivative of the self-energy to the vertex function (at least at one loop), and 
again show at the level of the Feynman integrands (with the proper momentum routing) that 


7NDR 
^mP 


(to = 0, A) = 1 . 


(3.10) 


Returning to the discussion of Z ]^ ttlce , we combine eqs. 
the chiral symmetry), giving 


and (pi), and use = 1 (since NDR preserves 


^lattice _ ^lattice (q 


(3.11) 


The corresponding relation holds for Zy tUce and Z^p lce . Note that all direct reference to the continuum has been 
eliminated. Since from now on we talk only about lattice quantities we can drop the qualifier “lattice” from eq. (3.11) 
and return to the not ation of Section [Hj. We remark that the renormalization factors were computed in Ref. ]lT| using 
the equivalent of eq. ( 3.11 ). 

The computation of the right hand side of eq. ( 3.1 1| ) is fairly easy, since it can be done in the massless limit. We 
need to compute the self-energy graphs (the continuunr-like graph and the lattice tadpole) and the vertex correction 
graph (for each current). To avoid computing any covariant parts, we route the external momentum through the 
fermion line. This routing will also be very convenient in the computation of q*. 

We write the lattice self-energy as 


S(p) = — E 0 (p 2 , to 2 , A 2 ) -H^E;l(p 2 ,to 2 , A 2 ) + mE 2 (p 2 , to 2 , A 2 ) . 
a 


(3.12) 


Only Ei is needed for wave function renormalization at to = 0; E 2 is needed for Z m p. Computing the self-energy 
graphs with the methods of Ref. Jh|] , we have 
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(3.13) 


£ 1 i2 (0,0,A 2 ) =g 2 C F [. 

Ji 


1,2 


Taking Wilson r = 1, the integrands I\, 2 are: 
1 1 


h = 


8Ai 4Ai A 2 

6»(tt 2 - £ 2 ) 


-2 + - 2 A 2 - + -J- (A 4 - As + A 1 A 4 + 4AiA 5 - 2 A 2 A 4 - 2 A 43 ) 

Z Z ZA2 


covariant part , 


4A 1 (A 1 -2) + A 4 A 4 — 2 4(9(tt 2 -£ 2 ) 

/! = si-+ 


2A 4 A 2 i A 
where A 4 ... A 5 are defined in Refs. p|,|l5| (for r = 1): 

A : =^sin 2 (y) , 

Ot 

A 4 = ^Sin 2 (4«) , 

a 


covariant part , 


(3.14) 


Ao — A4 “I - 4A? 


and we introduce 


Ai 3 = ^sin 4 (^A)sin 2 (£ Q ) . 


(3.15) 


(3.16) 


The “covariant parts” in eq. (3.14) do not need to be specified further, since they will cancel in the final answers, 
as explained above. The terms proportional to 1/f 4 , which come from the expansion in powers of a of added and 
subtracted covariant integrals, are written explicitly here, since they are needed to cancel the logarithmic divergences 
in the rest of the integrals. However, these terms will also cancel in the final results. 

In Ref. p| , the external momentum in the self-energy computation was chosen to flow through the gluon line. 
(This saved some effort in the expansion in powers of a, since the fermion propagator is more complicated than the 
gluon one.) For our current purposes, we therefore have performed the calculation of J 4 from scratch. In the case of 
/ 2 , however, one may set p = 0 from the beginning, so that the routing of the external momentum is irrelevant. Then 
I 2 may be taken immediately from Ref. [ p~ 4 | ]. 

For the vertex diagrams we may also set p = 0 ab initio and take the results from Ref. M]. We have 


where 


75,7m, 75 — (0V»75) 7/ii 75)5 Cf / 1^7^75 ,7m >75 

Jt 


>^7m75 ~~ 4 /C(, + 2 /Ci + /Co , 
Ji, = 4 /C 2 - 2 Kx + /Co , 
= I6/C2 - 4 /Ci + /C 0 


(3.17) 


(3.18) 


with 


IC' 2 
1C 1 
ic 0 


1 


A5 

48A| VAi 
1 “ 


- A. 


I6A1A2 

Y^j(a 4 + 4a l( A, -4)), 


9(tt 2 - e 2 ) 

4Z 3 


+ covariant part , 


12A 2 


A 4 — 4^ + 6 A 4 + 12A 2 

Ai 


(3.19) 


If we now write the Z factors in terms of the constants ca , cy and c m p as in eqs. ( |2. 13 2.1(j| ), we have 

ca = Cf J (Ii + C/7^75) , 
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c v =C F [{h + J llL ) , 

Ji 

c m p = C F [(h+J 75 ) ■ (3.20) 

Ji 


Explicit calculation shows that the covariant parts of I 12 an d 75j7(ii75 cancel in eg. ( 3.20|) , as expected from the 
general arguments above. Numerical integration of eq. (|3.20|) then reproduces eqs. (2.14, 


If). 


IV. TADPOLE IMPROVEMENT 


The idea behind tadpole improvement Jli] is based on the observation that, for the relatively large values of the 
lattice coupling constant g used in current lattice QCD computations, the expectation value of the link variable, 
uq = (U^ x ) (defined, e.g., in Landau gauge), is not very close to one. If one would calculate, for instance, the quark 
propagator in a mean-field approximation by replacing U^x —> uq in the quark action, one would find 

{q{x)q{y ))mf = Uo 1 (q(x)q(y)) {lee (M -> Uq 1 M) , (4.1) 

where the notation M —> Uq 1 M indicates that we replace M by u^ 1 M in the expression for the free Wilson propagator. 

This implies that, in order to make contact with the continuum, the quark field has to be renormalized by a factor 
y/uE. For local bilinears like A“(a;), P a (x) and V*(x), this results in a factor uq: 

Il^Ocont = ^0^^(*^)lattice j (4.2) 

and similarly for the other operators. From the mean-field result eq. @) for the quark propagator, we find that 
M c = Aruo and tomf = Uq 1 (M — M c ). (In terms of the hopping parameter n = 1/(2 M), we get n c = l/(8ruo).) 
Therefore, the operator (M — M c )P a (x) does not renormalize at all. 

The factor uq takes into account the gauge-field tadpole diagrams contributing to the renormalization constants Za 
and Zy. Of course, these are not the only contributions. We may now use standard perturbation theory to improve on 
our mean-field estimates for the renormalization constants, but we should be careful to take out those contributions 
that have already been absorbed into uq 0- It follows that 

^(a;)cont = -^ZaA^ x)lattice , (4.3) 

Uq 


where Uq T is the mean link ( U^ x ) calculated in perturbation theory, uq is the nonperturbatively ( i.e ., numerically) 
com pute d value (using the same definition!), and Za is the usual renormalization factor, calculated to one loop in 
eq. ([2.11: ). This leads to the definition of the tadpole-improved perturbative renormalization factor 


Za — Za/uq 


pt 


(4.4) 


with a similar definition of Zy. 

The mean-field link, u o, can be defined in different ways, and the idea of tadpole improvement is useful to the extent 
that the various definitions agree. Recently, Lepage |l7j and Trottier |L8]] have advocated the use of the Landau mean 
link as the nonperturbative uq factor. In Landau gauge, to one loop, we have (for N = 3) 



Uq T = 1 — 0.077466(l)g 2 (Landau gauge). 

(4.5) 

We may also calculate the 

mean link from M c : 



u p t = l- 0.108570(2) ff 2 (from M c ). 

(4.6) 

We may now write 

Zi = 1 — ag 2 , i = A, V , 

(4.7) 


and tabulate the numerical values for the constants Ci obtained from eqs. (4.4-4.7). 
gauge or from M c , we report the values in table 


For rtg T calculated in Landau 
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V. RESULTS FOR THE LEPAGE-MACKENZIE SCALE q 


The prescription for obtaining the Lepage-Mackenzie scale is to calculate q* from 

f d 4 qf(q) log {q 2 ) 


log((g*n = 


Jd 4 qf(q)) 


(5.1) 


if the one-loop integral for a particular renormalization constant has the form J d 4 qf(q ), and where q is the momentum 
flowing through the gluon line (see below). The condition that the momentum flow through the gluon line is critical 
both for the understanding and the correct application of the Lepage and Mackenzie scale-setting procedure. 

The Lepage and Mackenzie procedure is based on an older technique for setting the scale in continuum perturbative 
expressions due to Brodsky, Lepage, and Mackenzie (BLM) [|il|. In the BLM procedure, the scale is chosen to remove 
the leading-order coefficient of the number of fermion flavors from the perturbative expansion. Physically this absorbs 
a class of graphs from the fermion contribution to the vacuum polarization into the running coupling. This boils 
down to the requirement that q in eq. © is the momentum flowing through the gluon propagator for self-energy 
and vertex-correction diagrams. 

The original BLM scale setting procedure is more general than the Lepage and Mackenzie prescription (eq. ©)> 
because it can cope with perturbative expressions of operators with nonzero anomalous dimensions. The connection 
between the Lepage and Mackenzie scale setting procedure and the BLM procedure has been discussed by Kronfelcl Q. 

We will now calculate q*, both with and without tadpole impr ovement, using eq. ( |5.l[) . Let us, for definiteness , 
consider the case of Za- The relevant integral is given in eq. (3.2C), with I\ and 75 given by eqs. ( [3.14 3.IS,3.19). 
Since the loop momentum q has been chosen from the beginning to be the gluon momentum (Le., the external 
momentum is routed through the fermion line), we just need to recompute eq. ( {3.2C| ) with log(g 2 ) inserted. 

In table |J, we give the values of q* for the three Z{. and for the three cases of no improvement, Landau-gauge 
tadpole improvement, and M c tadpole improvement. For the case of Z m p, we did not include q* values from tadpole 
improvement, since the operator mP A (x) is not affected by it, as explained in the previous section. 

Note that if we had chosen the loop momentum to be the fermion momentum, the insertion of log(q 2 ) into the 
integrands would in general give a different result. While we may always shift the loop momentum in the original 
expressions for the diagrams, this is not true after we have inserted log ( q 2 ). Indeed, in the case where q is the fermion 
momentum, the integral with log(q 2 ) inserted is not even convergent. The covariant parts in that case generate 
integrals roughly of the type 


d 4 q 


2 2 

a z p z 


(q 2 + a 2 p 2 ) 3 


(5.2) 


which approaches a constant as a —> 0 but diverges logarithmically if log ( q 2 ) is inserted into the integrand before 
a —> 0 is taken. This tells us that this procedure of determining the typical scale q * is consistent at one loop only 
when this scale is the typical momentum of the gluon propagator. This is in accordance with the intuitive arguments 
underlying the BLM scheme. 

In previous computations of q* in the static-light case ^T| and in the Wilson or clover case ||. the loop momentum 
was not chosen to be the gluon momentum. The straightforward insertion of log( g 2 ) into the integrals was therefore 
not possible, because the “constant terms” (the a —> 0 limit of integrals like eq. (T2)) would give divergent results. 
These terms then had to be treated in an ad-hoc manner: Hernandez and Hill |21| replaced the integrals with constants 
over the Brillouin zone; Crisafulli et al. shifted and redefined the integration variables until an integral was found 
which remained convergent after insertion of log(g 2 ). Luckily, however, these terms contribute only a small amount 
to the final results for q*. Therefore the ambiguities in previous computations make only a small numerical difference. 
The results in Ref. j|] for the q* of Za or Zy differ by about 1% from the current, unambiguous answers. In making 
this comparison, we computed the integrals in Ref. |Q for ourselves (after correcting two typographical errors in I-p, 
in eq. (139)), since the quoted answers (q* = 2.4 for Zy and q * = 2.6 for Za without tadpole improvement) had too 
few significant figures to see the difference clea rly. 

Since the integrals which define the cq (eqs. ( |2 . 1 2 | , 2 . 16 ) or ( [ 3 . 20 | )) are finite and represent cutoff effects, one would 
anticipate q* for these quantities to be of order n/a. In the tadpole-improved case, where one hopes the most severe 
cutoff effects are taken out of the perturbative renormalization parts, one might expect the values of q* to be much 
lower Q]. Table FE] shows that, while q* is indeed reduced for tadpole-improved quantities, this reduction is rather 
minimal, and tadpole-improved q*' s are still substantially larger than 1/a. 





















VI. SCALE SETTING USING A NONPERTURBATIVE EXPRESSION FOR THE COUPLING 


The prescription for computing q* given in eq. (0 can be motivated starting from the assumption that a natural 
value for q* would be obtained from 


*( 9 *) J d 4 qf{q) = J d i qf{q)c 


'(<?) 


( 6 . 1 ) 


The idea here is to use a coupling that depends on the internal momentum in the diagram to weight the integr and , 
rather than multiplying the diagram with a coupling at a constant scale. Unfortunately the right hand side of eq. (0 
is singular if the perturbative renormalization group improved evolution equation is used for the coupling, because of 
the Landau singularity. Lepage and Mackenzie avoided this problem, by using the one loop perturbative evolution 
equation for the coupling (which is free of the Landau sing ularity) in eq. (0) to derive eq. (|5.l|). 

As a consistency che ck o n the scale obtained from eq. ( |5.l| ), it is instructive to try to use a nonperturbative definition 
of the coupling in eq. ( |p.l| ) . The Landau singularity in the renormalization group evolution equation for the coupling 
occurs when the coupling is large, and the perturbative derivation of the evolution equation breaks down. Neubert |2^ j 
developed a formalism based on eq. (0 to calculate the scale for a number of continuum perturbative expressions. 
He found that his more general formalism produced well behaved perturbative expression for some quantitie s th at 
were previously thought to be unreliab le b ecause of their low BLM scale (the continuum analog of using eq. (5.1)). 
This motivates studying the use of eq. (0 for lattice perturbation theory. 

As the running coupling of QCD has been obtained nonperturbatively, from a number of lattice QCD simula¬ 
tions J23| |25|, it seems natural to try a coupling measured from a lattice QCD simulation in eq. ( |6.1| ) to estimate the 
scale for lattice perturbative expressions. A useful expression for a nonperturbative QCD coupling for scale setting 
was suggested by Klassen |2(J. This QCD coupling is defined by fitting the static potential between two quarks in 
momentum space to an ansatz for the evolution of the coupling t hat has no Landau pole and reduces to the pertur¬ 
bative expressions in the weak coupling limit. Klassen used eq. (0 with his nonperturbative coupling to estimate 
the scale for the plaquette, which was consistent with the result from Lepage and Mackenzie’s formalism. 

Unfortunately, combining Klassen’s nonperturbative definition of the coupling with the integrands for Za in eq. (0) 
produces a divergent integral. The problem is that the integrands for Za behave like O(-i^) for small momentum 
q. Klassen’s coupling also behaves like O(^), thus producing a divergent result. The physical reason for the low 
momentum limit of Klassen’s coupling is the assumption that the static quark potential at large distances is linear in 
the quark separation. However, it is expected that at a certain separation between the quarks, string breaking will 
occur, which will cause the potential to flatten out. The crossover from a linearly rising potential to a flat potential 
has never been seen in lattice gauge theory simulations (see Refs. |2^^] for recent discussions). ^From lattice gauge 
theory, it is thus currently impossible to extract a coupling in momentum space, for use in eq. ( |6.l|) , that has the 
correct low momentum limit. 

There is some evidence from experiment that the low energy limit of the coupling tends to a finite limit in the small 
momentum limit. In the theoretical analysis of e + e _ event shapes, Dokshitzer and Webber j2tj] introduced the low 
momentum mean of the coupling. Their fit to experimental data yielded 


1 


r 2 GeV 


2 GeV . 


a(fM 2 )d/j, = 0.52 ± 0.04 . 


( 6 . 2 ) 


If a low momentum coupling with the limit of eq. (^(^) exists, then it would have a finite limit with the integrands 
for Za- 

Inspired by the work that obtained eq. ((12), Shirkov and Solovtsov f30|| derived an expression for the coupling 
without a Landau pole: 


1 


OZss [k ) — 


A 2 


/3b \log{k 2 /A 2 ) A 2 — k 2 


(6.3) 


Equation (6.3) is based on a dispersive theory and the one-loop RG evolution equation. The choice of scheme is 
reflected in the A parameter. This definitio n of the coupling, as well as its connection to renormalons, is reviewed in 
Refs. [ ^Ip2| . The analytic coupling in eq. (13) has been used to study the scale dependence in the continuum jh|. 

Although eq. (13) could just be substituted into the right-hand side of eq. & to provide an answer for the 
one-loop perturbative contribution, it is convenient to calculate a q* to compare with the scale obtained from the 
Lepage and Mackenzie scale. We define a q* from: 
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1 1 ^ /^/(g) (i/log(f^) + 

log((g*) 2 /A 2 ) + 1 - (<Z*) 2 /A 2 = fd*qf(q) ’ (6 ’ 4) 

where q and \ = aA are dimensionless. Using A = 0.14 (A = 0.0587), taken from Ref. ^(| at (3 = 6.0 (/? = 6.8), we 
obtain q* = 2.32 ( q* = 2.33) for Za , and q* = 2.10 ( q * = 2.10) for the tadpole improved (using M c ) Za • The two 
scales are both close to ones obtained from the Lepage and Mackenzie scheme. Again we see that tadpole improving 
the perturbative factor only causes a small reduction in the scale. Although eq. (^J) is less convenient to use with 
simulation data, because of the dependence on the scale through A, this formalism does provide a check on the Lepage 
and Mackenzie scale-setting procedure. 


VII. VISUALIZATION OF THE INTEGRANDS 


To try to understand why the estimate of the scale for tadpole improved perturbative factors was not close to the 
naive expectation -j, and to investigate the effects of having different weighting functions, we wanted to “visualize” 
the lattice integrands. Consider a lattice approximation to a continuum integral 



L — lL- 1 L — lL -1 


L 1 L -Li -J -L L 1 

Ja' 52'52'52'52 f( 2nx / L > 27r 2// I o 27T3/L, 27rf/ L) 

rc=0 y= 0 z—0 t= 0 

(7.1) 

jj^2f(QE)n(q E ) , 

(7.2) 


QE 


where qE is a single member of the equivalence class under hypercubic transformations. The function n(gs) gives the 
number of momenta in each equivalence class. This technique is the basis of an efficient method to integrate lattice 
Feynman d iag rams |34| . 

In figure |l| we plot the f(qE ) function for the Za integrand with and without tadpole improvement. The mag¬ 
nitude of the tadpole-improved integrand is much reduced over the original integrand. However, if we choose the 
normalization of the tadpole-improved integrand to equal the not improved integrand at a specific momentum, then 
the two figures agree qualitatively. Because the Lepage-Mackenzie scale does not depend on the overall normalization 
of the integrands, it is perhaps not surprising that the scales for the tadpole and nontadpole Za are so close. 


VIII. CONCLUSION 

Our main results for the perturbative scales are contained in table |fi|. Tadpole improvement does not significantly 
reduce q* for the operators considered here. In particular, the scales for the tadpole-improved operators are higher 
than the scale q* = 1/a suggested by the intuitive idea that tadpole improvement reduces the effect of the momentum 
near the cut off Jlj]. 

The closeness of the scale for the tadpole-improved Za to the scale for the tadpole-improved Za in the static 
limit ^T| ( aq* = 2.18) seems tp indicate that the mass dependence of q * is weak. This suggests that using a mass- 
independent q * with the results of the mass-dependent perturbative calculation of Za by Kuramashi |35f| is a good 
approximation. This was the procedure adopted in the calculation of fs by the MILC collaboration 
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no improvement 

Landau gauge 

through M c 

CA 

0.133373(2) 

0.055907(1) 

0.024803(1) 

C V 

0.174083(3) 

0.096617(2) 

0.065512(1) 

CmP 

0.081419(2) 

- 

- 


TABLE I. One-loop constants for tadpole-improved renormalization constants. 


quantity 

no improvement 

Landau gauge 

through M c 

Za 

2.533 

2.235 

2.316 

Zy 

2.370 

2.090 

2.051 

ZmP 

1.905 




TABLE II. The scale q * (in units of 1/a) for various quantities and improvement schemes. The errors are less than 1 in the 
last place. 
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FIG. 1. Plots of the integrand for Za with (diamonds) and without (crosses) tadpole improvement (through M c ). The 
squares are the tadpole improved integrand normalized to agree with the nontadpole improved numbers at a specific point. 
The lattice volume was 32 4 . For clarity only momentum with magnitude less than two are displayed. 




